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Abstract 



Collective behavior of pulse-coupled oscillators has been investigated widely. As an ex- 
ample of pulse-coupled networks, fireflies display many kinds of flashing patterns. Mirollo 
and Strogatz (1990) proposed a pulse-coupled oscillator model to explain the synchroniza- 
tion of South East Asian fireflies {Pteroptyx malaccae). However, transmission delays were 
not considered in their model. In fact, the presence of transmission delays can lead to 
desychronization. In this paper, pulse-coupled oscillator networks with delayed excitatory 
coupling are studied. Our main result is that under reasonable assumptions, pulse-coupled 
oscillator networks with delayed excitatory coupling can not achieve complete synchroniza- 
tion, which can explain why another species of fireflies {Photinus pyralis) rarely synchronizes 
flashing. Finally, two numerical simulations are given. In the first simulation, we illustrate 
that even if all the initial phases are very close to each other, there could still be big 
variations in the times to process the pulses in the pipeline. It implies that asymptotical 
synchronization typically also cannot be achieved. In the second simulation, we exhibit a 
phenomenon of clustering synchronization. 
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1 Introduction 



Fireflies provide one of tlie most spectacular examples of synchronization in nature [H El 
|3l H]. Through observation, people discovered that there are several synchronization patterns 
in different species of fireflies. Some species, such as Pteroptyx malaccae, flash rhythmically 
in perfect unison [1]. Some other species, such as Photinus pyralis, usually show clustering 
synchrony or wave sweeping synchrony, instead of complete synchrony [2] . 

Many kinds of biological models have been studied for flashing behavior of fireflies. A pioneer- 
ing work for synchronization of fireflies has been reported by Buck in [2] . He investigated a variety 
of firefly species including Pteroptyx malaccae, Pteroptyx cribellata and Photinus pyralis. He sug- 
gested two kinds of flash models, phase-advance entrainment model and phase-delay entrainment 
model. In the phase advance model displayed in Fig. 1(a), an excitation level of pace-maker in 
the firefly brain is enhanced by an external light stimulation. Each time excitation reaches a 
threshold level, a neural signal is transmitted into the light organ in the abdomen and then a 
flashing light is produced. An example of this model is Photinus pyralis, the whole group of the 
species rarely synchronizes flashing [2]. Instead, "wave", "chain" or "sweeping" synchrony has 
been reported in some species with this phase advance model. A cluster of male fireflies flash 
and then the flash of light is dispersed into neighbor males. In the phase delay model displayed 
in Fig. 1(b), an excitation potential of pace-maker is reset to the basal level by a light stimula- 
tion and the potential increase is restarted to reach the threshold level for flashing. Pteroptyx 
cribellata is known as an example of the phase delay model and it usually shows a complete 
synchroniz at ion . 
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Fig.l: (a) Phase advance model (b) Phase delay model. 
(Reprinted from ^) 



Fig. 2: Graph of the function /. The time-course 
of the oscillator is given by a; = /(</?)• 
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MiroUo and Strogatz (1990) argued that Buck's model is not appropriate for synchronization 
because he assumed a hnearly increasing potential toward threshold as shown in Fig.l. Inspired 
by Peskin's model for self-synchronization of the cardiac pacemaker [5], MiroUo and Strogatz 
proposed a pulse-coupled oscillator model with undelayed excitatory coupling to explain the 
synchronization of huge congregations of South East Asian fireflies [Pteroptyx malaccae) [6J. This 
model is a network of pulse-coupled oscillators. Each oscillator is characterized by a state 
variable Xi which is assumed to increase toward a threshold at Xj = 1 according to Xi = fifi), 
where / : [0, 1] [0, 1] is smooth, monotonic increasing, and concave down, i.e., /' > and 
/" < 0. Here (fi G [0, 1] is a phase variable such that (i) dipi/dt = 1/T, where T is the cycle 
period, (ii) = when the ith oscillator is at its lowest state Xi = 0, and (iii) yjj = 1 at the end 
of the cycle when the ith oscillator reaches the threshold Xi = 1. Therefore / satisfies /(O) = 0, 
/(I) = 1. Fig.2 shows the graph of a typical /. When Xi reaches the threshold, the ith oscillator 
"fires" and Xj jumps back instantly to zero, after which the cycle repeats. The oscillators are 
assumed to interact by a simple form of pulse-coupling: when a given oscillator fires, it pulls all 
the other oscillators up by an amount e, or pulls them up to firing, whichever is less. That is, 

Xi{t) = 1 ^ Xj{t^) = min(l, Xj{t) + e), Vj i. (1) 

The main result in [6] is that for all and for almost all initial conditions, the system eventually 
becomes completely synchronized. Here, the concept of complete synchronization is defined as 
follows: if there exists a to > such that 

Xi{t) = Xj{t), for all t >to and all i ^ j (2) 

then the pulse-coupled oscillator network is said to be completely synchronized or, for simplicity, 
synchronized. 

Since MiroUo and Strogatz's model was introduced, many results on pulse-coupled networks 
with undelayed coupling, including undelayed excitatory coupling and undelayed inhibitory cou- 
phng, have been obtained [3 [HI [9l [TOl [TTl [12]. However, reaction and transmission delays are 
unavoidable in real biological systems. For example, normally the transmission delays of most 
fireflies from sensors to motor actions of flashing are around 200ms (see [2]). In fact, transmission 
delays influence the performance of synchronization for both excitatory and inhibitory coupling 
[T3t [HI [15] . In particular, it has been shown that for excitatory coupling, the presence of trans- 
mission delays can lead to desychronization [131 HH UHl [II]- To the best of our knowledge, this 
desynchronization was proved only for the case of two pulse-coupled oscillators, while for the 
case of A^ > 2 pulse-coupled osciUators, it was revealed only in simulations. 
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In this paper, pulse-coupled oscillator networks with delayed excitatory coupling are studied. 
We retain two of Mirollo and Strogatz's assumptions: the oscillators have identical dynamics, 
and all oscillators are coupled to all the others. Our main contribution is that we prove that 
N > 2 pulse-coupled oscillators with delayed excitatory coupling can not achieve complete syn- 
chronization. This result can explain why Photinus pyralis rarely synchronizes flashing, which 
is known as an example of pulse-coupled oscillator networks with delayed excitatory coupling 
[2]. Finally, two numerical simulations are given. In Simulation 1, we illustrate that even if all 
the initial phases are very close to each other, there could still be big variations in the times 
to process the pulses in the pipeline. It implies that asymptotical synchronization typically also 
cannot be achieved. In Simulation 2, we exhibit a phenomenon of clustering synchronization. 

The rest of the paper is organized as follows: In Section 2, we describe the network model. 
In Section 3, some deflnitions and notations are given. In Section 4, some lemmas are proved. 
In Section 5, we prove the main theorem. Two numerical simulations are given in Section 6. We 
conclude the paper in Section 7. 



The network consists of > 2 pulse-coupled oscillators with delayed excitatory coupling. As 
in [6], the oscillators have identical dynamics, and all oscillators are coupled to all the others. 
Each oscillator is characterized by a state variable Xi which is assumed to increase toward a 
threshold at Xj = 1 according to = f{^Pi), where / : [0,1] —>■ [0,1] is smooth, monotonic 
increasing, and concave down, i.e., /' > and /" < 0. Here (pi G [0, 1] is a phase variable such 
that (i) dipi/dt = 1/T, where T is the cycle period, (ii) (pi = when the zth oscillator is at 
its lowest state Xi = 0, and (iii) yjj = 1 at the end of the cycle when the ith oscillator reaches 
the threshold Xi = 1. When Xi reaches the threshold, the ith oscillator flres and Xi jumps back 
instantly to zero, after which the cycle repeats. That is. 



The new features are that, (i) in order to make the discussion a little easier, we let T = 1, 
i.e., d(fi/dt = 1, and (ii) the transmission delay r is introduced into the model. Because of the 
presence of the transmission delay, the oscillators interact by a new form: when a given oscillator 
fires at time t, it emits a spike; after a transmission delay r, the spike reaches all the other 
oscillators at time t + r and pulls them up by an amount e, or pulls them up to firing, whichever 
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is less. That is, 

x,{t) = l^Xj{t + T)=mm{l,xj{{t + T)-)+e), \/j ^ i, (4) 

instead of ([I]). Of course, if there are m oscillators firing simultaneously, then the state variables 
X of the m oscillators are increased by an amount (m — l)e and the state variables x of all the 
other oscillators are increased by an amount me. By ([3]) and (jl]), it is clear that if a spike reaches 
the ith oscillator just at the moment that the ith oscillator reaches the threshold Xi = 1, then 
the state variable Xi changes as follows: 

Xi{t) = min(l, Xi{t~) + e) = min(l, 1 + e) = 1 ^ Xi{t^) = 0. 

In addition, we make the following two assumptions: 
(Al) The system is started at time t = with a set of initial states < Xi(0) < 1, and there are 

no firings in time [— r, 0). 
(A2) The transmission delay r and the coupling strength e satisfy f{2T) + Ne < 1. 

The assumption (Al) about the initial conditions is quite common for pulse-coupled networks 
with transmission delays [131 El HQ IE] • The time t = can be regarded as the moment at 
which fireflies assemble and begin to flash. Since jumps back instantly to when it reaches 
the threshold 1, the state values and 1 may be considered to be the same state. Therefore, we 
assume that the initial states satisfy < Xj(0) < 1 instead of < Xj(0) < 1. 

The assumption (A2) means that both r and e are relatively small. It is reasonable to fireflies 
with the phase advance entrainment model, especially to Photinus pyralis. The transmission 
delays of most fireflies from sensors to motor actions of flashing are around 200ms, while the 
regular endogenous period of Photinus pyralis is almost 6000ms (see |2j). Normally, excitatory 
activations generate frequent spikes in a short period of time, because it can shorten the period 
of oscillator by increasing the potential level (see [H]). It means that strong excitations result 
in frequent flashing of light, which will exhaust out the energy of the insect. Thus, a relatively 
small coupling strength for excitatory action is used in many species of fireflies, which follow the 
phase advance model. In the following sections, one can see that under this assumption, the 
number of case distinctions will be reduced. 

With the monotonicity of /, the state variable Xi and the phase variable (pi are one-to-one 
correspondence. Therefore, the synchronization of the state variables xi, . . . ,xn is equivalent to 
the synchronization of the phase variables (pi, . . . , ipi^i. In the following, instead of investigating 
Xi, we investigate dynamical behaviors of ipi directly. 
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Prom the above description, the phase variable processes the following properties. 



Proposition 1 The phase variable (pi satisfies: 

(a) Ifm spikes reach the ith oscillator at time t, then (pi{t) — /~^(min[l, +me]); 
If no spikes reach the ith oscillator at time t, then (fii{t) — (fii{t~). 

(b) Ifipiit) = 1, then ip,{t+)=0; 
Ififiit) < 1, then ifi,{t+) = ifi{t). 

(c) If no spikes reach the ith oscillator in time (^1,^2) and the ith oscillator do not fire in time 
{ti,t2), then ipiit^) = (Pi{tt) + {t2 -ti). 

Remcirk 1 The assumption (Al) and Proposition 1(b) imply that ipi(t) 7^ for all t > and 
all 1 < i < N . Thus, in our model, each phase variable (pi satisfies (fi{t) e (0, 1] for allt>0. 

3 Definitions and notations 

In this paper, the concept of complete synchronization is defined as follows: 

Definition 1 // there exists a to > such that (fi{t) — (fij{t) for all t > to, then we say that 
oscillator i and j can achieve complete synchronization or, for simplicity, synchronization; and 
say that oscillator i and j have already been completely synchronized at time to or, for simplicity, 
synchronized. 

Definition 2 // there exists a to > such that ^pi{t) = ^2{t) = ■ ■ ■ = <fN{t) for all t > to, 
then we say that the pulse-coupled oscillator network can achieve complete synchronization or, 
for simplicity, synchronization; and say that the pulse-coupled oscillator network has already been 
completely synchronized at time to or, for simplicity, synchronized. 

Throughout the paper, the following notations will be used. 

(i) FmiO) = f-\mm[l, f{d) + me]), where m e Z+ and < ^ < 1. Here, Z+ = e > 0}. 

(ii) Tj = {t G M+|t is the time at which the ith oscillator fires}, i = 1, . . . , A^. Here, = {2; e 
R|z > 0}. 

(iii) Diit) ^{q^t' + T -t\t' eTi^^[t-T, t)}, i ^ 1, . . . , N, t > 0. 

(iv) Di{t+) ^ {rj ^ t' + r - t \ t' e Ti n {t - T,t\}, i ^ 1, . . . , N, t > 0, which is called the 
right-limit of Di{t) at time t. 
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Remark 2 According to Proposition 1, if m spikes reach the ith oscillator at time t, then 
(Pi{t) — /~^(min[l, +me]) — Fjn{(fi{t~)); if no spikes reach the ith oscillator at time 
t, then (fii{t) — (piit") — Fo{(fi{t~)) . Therefore, for the convenience of later use, we introduced 
the notation Fm{0). 

Remark 3 The notation Ti is called the firing-time set of the ith oscillator. Using the firing- 
time set Ti, we defined two other notations Di{t) and Di{t^). They can be used to calculate the 
reaching time of the spikes which the ith oscillator has emitted, but the other oscillators have 
not received. From the definition of Di{t), one can see that for each rj e Di{t), there must be a 
spike of the ith oscillator, which will reach the other oscillators at time t -\- rj. In particular, if 
e Di{t), then there must he a spike of the ith oscillator, which reaches the other oscillators at 
time t. 



Remark 4 The notations Di(t) and Di(t~^) will play an important role in our discussion. For 
the pulse-coupled networks without transmission delays, instantaneous synchronization (v5i(to) = 
ifj{tQ)) shows that oscillator i and j have been completely synchronized at time to {^Pi{t) = <fj{t) 
for allt > to), because the coupling is all-to-all and the oscillators have identical dynamics. But, 
for those with transmission delays, this criterion fails {instantaneous synchronization does not 
mean complete synchronization) , because of the presence of delays. Fig. 3 shows an example in 
which the criterion fails. One key role of Di{t) and Di{t'^) in this paper is to judge whether the 
pulse-coupled oscillator network with delays has been completely synchronized at time to, when 
we have >^i{to) — . . . — 9?Ar(to)- 
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(p T — 4> T T + <p 

(d) (e) (f) 

Fig. 3: A system of two oscillators (A and B) governed by a; = /(i^), with r and e satisfying e < /{t) < 1. 
The initial phases are chosen to be <^a(0) = 1 — and <^s(0) = 1 such that < <p < t and /(r — (p) + s = /(r). 
(a) The state of the system at t = 0. At the moment oscillator B emits a spike, (b) The state of the system just 
after t = 0. (c) The state of the system a.t t = (p. At the moment oscillator A emits a spike, (d) The state of the 
system just after t = cf). (e) The state of the system at t = t. At the moment the spike emitted by B at f = 
reaches A. (f) The state of the system at t = t + (p. At the moment the spike emitted hy A at t = (f) reaches B. 
Although equality <^a(s) = fsis) holds for all s G [t,t + 4>), (f) shows that the system has not been synchronized 
at t = T + <j). 

Proposition 2 Fm{0) has the following properties: 

(a) < F^{e) <1 for alio <e <1 and all m e Z+. 

(b) If 9i < 02, then Fm{Oi) < Fm(^2)- In particular, if 9i < 62 and Fjn{Oi) < 1, then 
Fmidi) < Fm{02)- {the monotonicity with respect to 6) 

(c) //mi < m2, then Fm^{6) < F^^^O). In particular, if mi < m2 and F^^iO) < 1, then 
Fmi{0) < Fm2{d)- {the monotonicity with respect to m) 

(d) When < Fm{0) < 1, the derivative of Fm{0) with respect to 9 exists, and ^^^^^^ > 1- 

(e) Ife,5>0 and Fm{e + S)<1, then F^(^) + 5< F^{e + 5) 

and the equal sign holds if and only if 5 — or m — 0, i.e., m5 — 0. 

(f) Fn{Fm{9)) = Fm+n{0) for all < 9 < 1 and all m,neZ+. 

(g) If 9, . . . , 5„ > and F^,+...+^„(6' + 5i H h (5„) < 1, then 

Fm„ {■ ■ ■ (F„,2 {F„i^ {6) + 61) + 62) -\ ) + Sn< Ffni-i h'Tin 

{d + 6i + --- + 

and the equal sign holds if and only if mid i = ■ ■ ■ = mnSn = 0. 

Remcirk 5 Property (d) comes from the concavity assumption of function f. Property (e) 
follows from {d), (f) is clear, and (g) follows from (e) and (f). 
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Propositions Di{t) and Di{t~^) have the following properties: 

(a) For each rj e Di{t), we have <r) < t. 
(a') For each rj e we have < 77 < r. 

(b) For any r]i,r]2 & Di{t), we have \r]i — 772] < t. 
(b') For any rji,rj2 G Di(t^), we have \r]i — 7/2I < t. 

(c) Ifrje Diit), then t - [t - rj] e Ti, i.e., (fi{t - [t - 7]]) = 1. 
(c') Ifrje Di{t+), then i - [r - 77] e T^, i.e., ipi{t - [r - rj]) = 1. 

(d) IfrjE Di{t) and < rj <r, then rj e Di{t^) too. 
(d') IfrjE Di{t'^) and <r] <t, then rj e Di{t) too. 

(e) If f] E Diit), then rj — cu E Di{t + ui), where cu G [rj — T,rj]. 

(e') If T} E Di{t~^), then rj — uj E Di{{t + co")"*"), where uj E [rj — T,r)). 

4 Properties of the pulse-coupled oscillator network with 
delayed excitatory coupling under the assumptions (Al) 
and (A2) 

It is difficult to prove the main theorem directly. The following lemmas are needed. Lemma 
1 shows that under the assumptions (Al) and (A2), the time between any two firing activities 
of any oscillator is longer than 2r. This result is reasonable to fireflies, because the transmission 
delays of flreflies arc relatively short and flreflies can not flre twice in such a short period of time. 
By Lemma 1, we obtain Lemma 2, which shows that under the assumptions (Al) and (A2), both 
Di[t) and Di[t~^) have at most one element. In Lemmas 3-5, three criteria of synchronization of 
two oscillators are given. In Lemma 6 and Lemma 7, two special cases are discussed, respectively. 
The two lemmas show that the network can not achieve complete synchronization even if it has 
been divided into two cliques. 

Lemma 1 // an oscillator fires at time ti and t2 with ti ^ t2, then \ti — ^2! > 2t. 

Proof: Assume ^2 > and that the oscillator is the ith oscillator. Then we have ^Pi{ti ) = and 
^iih) — 1- For ti,t2, there must be ^3,^4 G [^1,^2] such that the ith oscillator fires at ^3,^4 and 
does not fire in (^3, ^4). If we can prove 1^3 — ^4] > 2t, then \ti — ^2! > 2t must hold. Therefore, 
without loss of generality, we can assume that the ith oscillator does not fire in time (ti, ^2)- 
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Suppose t2 — ti < 2t. 
Case 1: t2 < r. 

Because of the presence of the transmission delay and the assumption (Al), no spikes can 
reach the ith oscillator in time (^1,^2] Q (0,t). Then, according to Proposition 1, (pi{t2) — 
^1(^2) — Viiti) + (^2 — ti) — t2 — ti < T < 1, which contradicts that (pi{t2) — 1. So, this case is 
impossible. 

Case 2: t2 = t. 

Because of the presence of the transmission delay and the assumption (Al), only the spikes 
that the other oscillators emit at time can reach the ith oscillator in time (ti, t2]- Assume that 
among the other oscillators, there are m oscillators firing at time 0. (Here, m can be equal to 
zero. It means that among the other oscillators there arc no oscillators firing at time 0. For the 
case of m = 0, the following proof holds true.) It is obvious that < m < — 1. According to 
Proposition 1 and the assumption (A2), we have 

<^i(^2 ) = ^iitt) + {h - h) ^ T - h < T 

and 

'Pi{t2) = f {inm[l, fit^iit^)) + me]) 

= f-\mm[lj{t2-h) + me]) 

< r\mm[l,f(r) + (N-l)e]) 

< r\mm[lJ{2T) + Ne]) 
= r\f{^T) + Ns) 

< /-'(!) = 1 

which contradicts that (pi{t2) = 1. So, this case is impossible. 
Case 3: ^2 > 'T- 

Assume that there are m = mi + m2 + • • • + m; spikes which reach the ith oscillator in time 
{ti,t2\: 

mi spikes reach at time ti + 5i; 
m2 spikes reach at time h + 61 + 62; 



mi spikes reach at time ti + Si + ■ ■ • + Si. 



9 



where Si, . . . ,Si > 0, and Si + ■ ■ ■ + Si < t2 — ti. (Here, m can be equal to zero. It means that 
there are no spikes reaching the ith oscillator in time (^1,^2]- When m — 0, we can let I — 1, 
nil — and Si be any value in the range (0, t2 — h] . Hence the following proof holds true for the 
case of m = 0.) 

Because the ith oscillator docs not fire in time (^1,^2), according to Proposition 1 we have 
fiiih + Si]-) = (p,{t+) + Si =61 
^i{ti + Si) = FmMih + Si]-)) = FmdSi) 
y:>i{[ti + Si+ 52]-) = Viiti + 5i) + 52 = FmdSi) + S2 
Vi{ti +5i+ 82) = Fm,{<Pi{[ti 62]-)) = Fm,{Fm,iSi) + S2) 



+5i + --- + 5i]-) = ipiih + <5i + • • • + Si-i) + 5i = Fm,_, (■ • • Fm, {5i) + ^2) • • • + ^i-i) + 5i 
Viitl +5i + --- + 6l)=Fm, Mti +6i + --- + 5i]-)) = Fm, {Fm,_, (• • • Fm,iFm, (Si) + ^2) • • • + Sl-i) + 61) 

Because of the presence of the transmission delay and the assumption (Al), the spikes which 
reach the ith oscillator in time {ti,t2\ are the ones that the other oscillators emit in time {ti — 
'T, ^2 — t] n [0, +00) . If all the other oscillators fire not more than once in {ti — r, ^2 — t] fl [0, +00) , 
i.e., m < N — 1, then it follows from the assumption (A2) and Proposition 2(g) that 

^i{t2) = <Pi{tl +Si + --- + Si) + [t2-'{ti+Si + --- + Si)] 

= Fm, (F„,_, (• • • F^,{Fr„,{Si) + (52) • • • + + 61) + [t2-{ti+Si + --- + 61)] 

< Fm,+...+m, {6l + --- + Sl + [t2- {tl+Si + --- + St)]) 
= Fmi-\ \-mi {t2 — tl) 

< Fn{2t) < 1 

which contradicts that ^pi{t2) = 1. Hence there must be some oscillator firing more than 
once in {ti — T,t2 — t] (1 [0, +00). It means that there must be time ^3,^4 with ^3 < ^4 and 
^3, ^4 e {tl —T,t2 — T]r\ [0, +00) such that some oscillator fires at t:^, t^^ and does not fire in (t3, ti). 
If ^4 > T, as above, there must be time ^5, t^ with ^5 < t^ and ^5, t^ e (^3 — r, ^4 — r] fl [0, +00) 
such that some oscillator fires at ^5, t^ and does not fire in [t^, to). litG>r, ■ ■ • This leads to a 
finite sequence: 

(tl, ^2) — ^ (^3) ^4) — ^ ■ • • — ^ {t2k-l, t2k) {t2n-li t2n) 

which satisfies that 

t2(k-i)-i — T < t2k-i < t2k < ^2(^-1) — r. A; = 2, . . . , n, 

t2k-i > 0, A; = 1, . . . ,n, 

t2k > T,k ^ 1, . . . ,n - 1, and t2n < r. 
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The last term of the sequence means that some oscillator fires at t2„-i and and does not fire 
in (i2n-i,^2n), with < ^271-1 < ^2n < By case 1 and case 2, this is impossible. So, case 3 is 
impossible, too. 

Therefore, we have ^2 — > 2r, which completes the proof. □ 
By Lemma 1, we get Lemma 2. 

Lemma 2 Under the assumptions (Al) and (A2), Di{t) and Di{t^) have the following prop- 
erties: 

(a) Di{t) — {77} with <r] < T or Di{t) — $, for all 1 < i < N and all t > 0; i.e., set Di{t) has 
at most one element for all 1 < i < N and allt>0. 

(b) Di{t+) = {t]} with <r]<T or Di{t+) = 0, for all 1 < t < N and all t > 0; i.e., set Di{t+) 
has at most one element for all 1 < i < N and all t >0. 

(c) IfDi{t) ^ 0, then 0<ipi{t) <1. 

(d) IfDi{t) = {771}, Dj{t) = {772} and 771 < 772, then (pi{t-) > (pj{t~). 

Proof: (a) Assume that there are 771, 772 € Di{t) with 771 7^ 7^2- Then by Proposition 3(c) we have 
(fi{t— [t — 771]) — (fi{t — [t — r]2]) — 1, i.e., the ith oscillator fires at time t— [t — 771] and i — [r — 772]. 
However, it follows from Proposition 3(b) that |[^ — (t — 771)] — [t — {r — 772)]] = I771 — 772I < r. 
This contradicts Lemma 1. So, for any 771,772 G Di{t), the equality 771 = 772 must hold. It shows 
that Di{t) = {7;} with < 77 < r, or Di{t) = 0. 

(b) The proof is similar to that of (a). 

(c) If Di{t) 7^ 0, there must be a 7; e A(i) with < 7; < r. By Proposition 3(c), we have 
(fi{t— [t — rj]) = 1. Since t — [t — [r — rj)] = t — rj < 2r, by Lemma 1, the ith oscillator never 
fires at time t. It means (pi{t) < 1. Thus, from Remark 1, we have < (pi{t) < 1. 

(d) Assume that there are m — mi + m2 + ■ ■ • + mi spikes which reach the ith oscillator in 
time {t — {t — 772), t): 

mi spikes reach at time t — (r — 772) + 

77T.2 spikes reach at time t — (r — 772) + 5i + ^2; 

mi spikes reach at time t — {t — r]2) + Si + ■ ■ ■ + Si. 

where Si, . . . ,Si > 0, and Si + ■ ■ ■ + Si < r — 772. (Here, m can be equal to zero. It means that 
there are no spikes reaching the ith oscillator in time {t — {t — 772), i). When ttt. = 0, we can let 
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/ = 1, m; = and 5i be any value in the range (0, r — 772). Hence the following proof holds true 
for the case of m = 0.) 

Since Di{t) = {r/i}, Dj{f) = {772} and rji < 772, by Proposition 3(c) and 3(c) we have that 
(Pj{t - (r - 772)) = 1, Di{s) = {tji + {t-s)} for all s G [t - (r - 7j2),t] and Dj{s) = {7J2 + {t - s)} 
for all s e (t — (r — 772), i]. 

Prom Di{s) ^ {rji + {t- s)} {0} and Dj{s) ^ {r]2 + {t- s)} {0} for aR s e {t- {r -rj2),t), 
it follows that oscillator i and j can not receive the spikes of each other in time {t — {t — r)2),t), 
which means that the spikes reaching the jth oscillator in time {t — (r — 772), t) are the m spikes 
above. 

Because A(s) = {vi + - «)} 7^ and Dj{s) = {772 + {t-s)}y^$ for all s E (t - (t - r]2),t), 
from Lemma 2(c), it follows that < (pi{s) < 1 and < (pj{s) < 1 for all s e (t — (r — 772), i), 
that is, both oscillator i and j do not fire in time (t — (r — 7/2), t). 

Since Di{t — (r — 772)) = {t + 771 — 772} 7^ 0, by Lemma 2(c) we have < (pi{t — (r — 772)) < 1. 
Denote (p — (fi{t — (r — 772))- 

Therefore, according to Proposition 1, we have 

V>i{[t -{r- m) + <5i]") = 9^ + 

^Pi{t - (r - 772) + h) = Fm, {<i> + ^i) 

V>i{[t - (t - m) + 82]-) = Fm, {(f> + Si) + 62 

iPi{t - (r - 7?2) + Si + 62) = Fm,{Frm {<P + ^l) + ^2) 



ipi{[t - (t - r?2) + ^1 + • • • + 5i]-) = Fm,_, (• • • Fm,{Fm, {(j) + 5^) + 82) ■ ■ ■ + 5l-x) + 5l 

<fi{t-{T-Tj2) + Si + --- + Si) = Fmi{Fmi_,{- ' ' Fm2{FmA4> + Si) + S2) ■ ■ ■ + Si-i) + Si) 

and 

ip,i[t-{T-T^2)+Sir)^5l 

'Pji't - (t - m) + Si) = FmJSi) 

'Pj{[t-{T-V2) + Si+S2r) = F^,{Si)+S2 
<fj{t - {t - ri2) + Si + S2) = Fm^ {Fm, {Si) + S2) 



Vj{[t -iT-r]2)+Si + --- + Si]-) = Fm,_, (• • • Fm, {Fra, (^i) + ^2) • • • + Sl-i) + Si 
iPj{t -{T-r]2) + Si + --- + Sl) = Fm, {Fm,_, (• • • Fm, {Fm, {Si) +S2)--- + Sl-i) + Si) 
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By Proposition 2(b), we have 

ifiit - (r - + Si) = Fm^ {(/) + Si) > Fm, {Si) = ipj{t - (r - rj^) + Si) 

ifiiit - {t - r]2) + Si + S2) = Fm, (0 + ,5i) + S2) > Fm^ (F™, {5i) + S2) = (fij {t - {t - m) + Si + S2) 



<Pi{t - (r - %) + <5i + • • • + Si) = Fm, (• • • F^, {Fm, {(I) + Si) + S2) ■ ■ ■ + Si) 

> Fm,{- ■ ■ Fm,{FmA^l) + ^2) ■ ■ ■ + Si) = ifijit - (j - 1^2) + Si + ■ ■ ■ + Si) 

Therefore, 

ipi{t-) = ipi{t-{T-r]2) + 5i + --- + 5i) + {t-[t-{T-r]2) + h + --- + 5i\] 

> ipj{t -{r-r]2) + 6i + --- + 6i) + {t-[t-{T-r]2) + 6i + --- + 6i]} 

= ^jin □ 



In the following (Lemmas 3-5), three synchronization criteria of two oscillators will be given. 

Lemma 3 Oscillator i and j have been synchronized at time to, if and only if the two oscilla- 
tors satisfy (fii{to) = fjito) and A(4) = Dj{t^). 

Proof: (1) Suppose that ipi{to) = ^jito) and -Di(^o ) = Djit^). 

It means that the two oscillators have acted as one, because their dynamics are identical and 
they are coupled in the same way to all the other oscillators. That is, oscillator i and j have 
been synchronized at t^. 

(2) Suppose that oscillator i and j have been synchronized at namely, (pi{t) — (pj{t) for all 
t > to. 

We only need to prove A(^o') = -Dj(^o )• Assume A(^o ) 7^ A(^o )• 

It is clear that A(^o ) 7^ or Dj{t^) ^ 0. Without loss of generality, we let A(io ) 7^ 0- By 
Lemma 2(b), we have A(^o ) — {v} -^^j(^o) '^^^^ < 77 < r. Then by Proposition 3, we get 
A(^o + v) — {0} Dj{to + 77), which means that a spike of oscillator i will reach oscillator j at 
time to-\-r], but no spikes of oscillator j will reach oscillator i at time to + r]. We let be the 
number of the spikes which reach oscillator i at time to + and ruj be the number of the spikes 
which reach oscillator j at time + V- So, the equality rrii = rrij — 1 holds, because the two 
oscillators arc coupled in the same way to all the other oscillators. It is obvious that 99,; ([to +??]") = 
fj{[to + because (pi{t) — (pj{t) for all t > to- Then, according to Proposition 1, we have 
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Vi{to + v) = f ^(min[l,/((^j([to + r;] )) + mi€]) ^ f ^{mm[l, f{(pj {[to + 1]] ))+mje]) = ipj{to + ri), 
unless f{(fi{[to + r]]~))+mie > 1 and f{ipj{[to + r]]~))+mje > 1, i.e., (pi{to + r]) = ifij{to + r]) = 1. 
However, by Di{to+ri) — {0} 7^ and Lemma 2(c), we have (pi{to+ri) < 1. This is a contradiction. 
Thus, Di(4) = Dj(t+). □ 

Lemma 4 // Di{to) = Dj{to) — {77} with < 77 < r, t/^en oscillator i and j have been 
synchronized at time to — {t — rf); conversely, if there exists ati < to such that oscillator i and 
j have been synchronized at time ti, then Di{to) — Dj{to). 

Proof: (1) Suppose that A(^o) = Dj{to) = {?]} with Q <ri <t. 

Prom Proposition 3, we have Di{[to—{T—r])]'^) — Dj{[to—{r—r])]'^) — {r} and ipi{to—{T—r])) — 
V'ji'to ~v)) — 1- By Lemma 3, oscillator i and j have been synchronized at time to — {t — rj). 

(2) Suppose that there exists ati < to such that oscillator i and j have been synchronized at 
time ti. 

Assume A(^o) 7^ Dj{to). It is clear that Di{to) 7^ or Dj{to) 7^ 0. Without loss of generahty, 
we let -Dj(to) 7^ 0- By Lemma 2(a), wc have Di{to) = {f]} 7^ Dj{to) with < rj < t. If follows 
from Proposition 3 that Di{to + i]) = {0} 7^ Dj{to + rj). Since oscillator i and j have been 
synchronized at time t\ with ti < to, by an argument similar to that used in the second part of 
the proof of Lemma 3, we can obtain (fi{to+r]) — (fj{to+r]) — 1. However, by Di{tQ+r]) — {0} 7^ 
and Lemma 2(c), we have (pi{to + rj) < 1. This is a contradiction. Thus, A(io) — Dj{to). □ 

Lemma 5 Oscillator i and j have been synchronized at time to > and they have not been 
synchronized at time t for all < t < to, if and only if the two oscillators satisfy one of the 
following cases: 

(a) (pi{to) 7^ <Pj{to), Di{to) = Dj{to) = and (pi{to) = (pj{to) = 1. 

(b) cpiito) 7^ ^j{to), Di{to) = {0}, Dj{to) = and ipi{to) = ipj{to) < 1. 

(c) (pi{to) 7^ ^j{to), Di{to) = 0, Dj{to) = {0} and ipi{to) = ^j{to) < 1. 

Proof: (1) Suppose that oscillator i and j have been synchronized at to and that they have not 
been synchronized at t for all < t < to- 

Wc have (pi{to) = <fj{to), because oscillator i and j have been synchronized at to- Assume 
that (pi{to) = <^i(^o )• From Lemma 1, there must be a 5 > such that no oscillators fire and no 
spikes reach in time {to — S,to)- It means that d(fi{t)/dt — d(pj{t)/dt — 1 for all t & {to — S,to)- 
Thus, we get ipi{t) — (pj{t) for all t >to — S/2. This contradicts that oscillator i and j have not 
been synchronized at t for all < t < to- So, (pi{to) 7^ (pj{to)- 
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If no spikes reach both oscillator i and j at time Iq, then it follows from v?i(to ) 7^ ^j{to) ^^^^ 
^iito) — fii'to) 7^ fjiio) — fji'to)- It contradicts that oscillator i and j have been synchronized 
at to. Therefore, there must be a spike which reaches at least one of the two oscillators at time to, 
i.e, there must be an oscillator k such that Dk{to) — {0}. (This is an important property of pulse- 
coupled networks.) Let X be a set of oscillators: oscillator k' & K \l and only if Dk'{to) — {0}. 
And let m be the number of the oscillators in K. 

Case 1: oscillator i,j ^ K. 

Then f-\mm[lj{<fi{t^)) + me]) = <fM = fjito) = r\mm[l, f{<f^{t^)) + me]). But 
(Piit^) ^ (PjitQ). So, (pi{to) = (Pj{to) = 1. By Lemma 2(c), we get A(^o) = Dj{to) = 0. Then 
case (a) holds. 

Case 2: oscillator i e K and oscillator j ^ K. 

Then we have D,{to) = {0} ^ Dj{to). A(^o) = {0} shows that D,(t^) = and v?i(to) < 1- 
Since oscillator i and j have been synchronized at to, it follows from Lemma 3 that -Dj(io ) — 
Di{t^) = 0, which means Dj{to) ^ $ or {0}. But Dj{to) ^ {0}. So, D^{to) = 0. Then case (b) 
holds. 

Case 3: oscillator i ^ K and oscillator j e K. 
It is similar to case 2. And case (c) holds. 
Case 4: oscillator i,j e K. 

Then Di{to) — Dj{to) — {0}. By Lemma 4, oscillator i and j have been synchronized at time 
to — T < to, which contradicts that oscillator i and j have not been synchronized at t for all 
<t < to. Therefore, this case is impossible. 

(2) Suppose that one of the three cases (a), (b) and (c) holds. Then we have ipi{to) — <Pj{to) 
and A(^o ) — ^ji^o)- F^om Lemma 3, oscillator i and j have been synchronized at to. If there 
exists a t' < to such that oscillator i and j have been synchronized at t', then ^iit^) — (fij{tQ). 
This is a contradiction. So, oscillator i and j have not been synchronized at t for all < i < to.n 

Two special cases are discussed in the following two lemmas, which show that the network 
can not achieve synchronization even if it has been divided into two cliques. 

Lemma 6 // there exists a to > such that (pi{to) — ■ ■ • — (fm{to) — 1 — (f{0<(f<l), 
Di{to) = • • • = D^{to) = and (Pm+i{to) = • • • = </?jv(io) = 1; D^+i{to) = • • • = -Djv(^o) = 0; 
where 1 <m < N , then the pulse- coupled oscillator network can not achieve synchronization. 
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Proof: Let n — N — m. 



V^i{to) = • • • = ^Pm{to) = 1 - (0 < < 1), Di{to) = • • • = Dm{to) = and iPm+i{to) = • • • = 
V^Af(io) = 1, Dm+i{to) = • • • = -DAr(to) = show that oscillator 1, . . . , m have been synchronized 
at to, and oscillator m + 1, . . . ,N have also been synchronized at to- Therefore, we only need to 
discuss two oscillators, the ith oscillator and the jth oscillator [1 < i < m, m + 1 < j < N). The 
synchronization of the whole network is equivalent to the synchronization of oscillator i and j. 

According to the model rules, we get the following tables of ipi{t), Di{t), ^j{t) and Dj(t): 
Case 1: cp < r. 



time t 


fi{t) 


D^{t) 




Dj{t) 


to 


1-ip 





1 





4 










{r} 


to + V 


1 







{T-ip} 


{to + <P)+ 





{r] 


V> 


{T-if} 


to T 




{r) 




{()} 


ito+r)+ 


Fn{T-ip) 


W} 


Fn-l{T) 





to+T + ip 


Fm-l{Fn{T - if) + if) 


{0} 


Fm{Fn-i{r) + ^) 





{to+T + ip) + 


Fm-l{F„{T - if) + if) 





Frn{Fr,-i{T) + ^) 





to+T + ip+ 

1-F„(F„_i(t)+(p) 


l-[i^m(K-l(T)+^) 
-Fm-l{Fr,{T - if) + if)] 





1 






Case 2: (p >t. 



time t 


(Pt{t) 




'Pj{t) 


Dj{t) 


to 


1-^ 





1 





t^ 


1-^ 








{r} 


to+T 


F„(1-</? + t) 





Fn^l{T) 


{0} 


{to + r)+ 


Fn{l-(p + T) 





F„-i(t) 





to+T+ 

l-K(l-V3 + r) 


1 





1 - [F„il - if + t) 
-Fn-i{r)] 






Let 



Gi{9) = F„(F„_i(t)+^)-F^_i(F„(t-^) + ^) 0<e<T 

G2{e) = Fn{Fm-i{T) + e)-Fn-i{F^{r-e) + e) o<e<T 

Gs{e) = Fn{l-e + r)-Fn-i{r) r<e<l 

G,{e) = F^(1-^ + t)-F^_i(t) t<^<1 
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Define map of tlie form 



{Gi{9),P,Q) if <e <T,P = m,Q = n 

{G2{e),P,Q) if 0<e <T,P = n,Q = m 

{Gs{9),Q,P) if T<9<l,P = m,Q = n 

(G'4(^),Q,P) if T <e <1,P ^n,Q 

where 0<9<1, P — morn, Q — morn. 

Then the difference between the phase variables (pi and (fj can be shown by the following 
sequence: 

{if, m, n) —>■ G{ip, m, n) —>■ G'^{ip, m,n) —>■ ■ ■ ■ 

If the pulse-coupled oscillator network can achieve synchronization, then there must be a p G N 
such that GP{ip,m,n) = (0, m,n) or (0, n, m). In the following, we will prove that such p does 
not exist. 

When < e < T, we have < F„(t -9) <1,0< F^(F„_i(t) + ^) < 1 and < F„_i(F„(t - 
9) +9) < 1. Prom Proposition 2(d), it follows that F^(r - ^) > 1, F^(F„_i(r) + 9) > 1 and 
Pln-AFnir -9) + 9)>l, for all < ^ < r. Therefore, 

- F'm{Fn-.{r) + 9)- F'^_,{F^{t - 9) + 9) ■ [-F^(r - ^) + 1] > 1 

for all < ^ < T. Then from G'i(O) = 0, it follows that Gi{9) > for all < ^ < r. 

When T <9 <1, since 1 — 9 + t > t and F„_i(r) < 1, by Proposition 2(b) and 2(c) we have 
that 

G^{9) = Fn{l -9 + t)- Fn-l{T) > 

for all r < e < 1. 

Similarly, we can get that G2{9) > for all < ^ < r and G4{9) > for all r < ^ < 1. 

Thus, there docs not exist p G N such that G'\(p,m,n) = (0,m, n) or (0,n, m). It means 
that the pulse-coupled oscillator network can not achieve synchronization. □ 

Lemma 7 // there exists a to > such that (pi{to) — ■■■ — (fmito) — ^ (0 < < 1); 
Di{to) = • • • = -Dm(io) = (0 < 77 < r) and (pm+i{to) = • • • = <PN{to) = 1, -Dm+i(^o) = • • • = 

-Div(^o) = 0; then the pulse-coupled oscillator network can not achieve synchronization. 
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Proof: Similarly as in the proof of Lemma 6, we only need to discuss two oscillator i and j, 
with 1 < i < m and m + 1 < j < N. 

Let n — N — m. 

By -Di(to) = • • • = A-i(^o) = A+i(^o) = • • • = Dm{h) = {77} (0 < 77 < r) and Lemma 
1, oscillator 1, . . . , i — 1, i + 1, . . . , m do not fire in time [max(0, — {t — rj) — t), max(0, — 
r)]. By </?m+i(^o) = • • • = fNito) — 1 and Lemma 1, oscillator m + 1, . . . ,N do not fire in 
[max(0, to — {t — 7]) — t), max(0, to — t)], too. Therefore, no spikes can reach the ith oscillator 
in time [to — (t — r]),to]. Then we have V'i(to) = t — r], i.e., rj = t — ip. 

Now we have the following table of (pi{t), Di{t), (pj{t) and Dj{t): 



time t 


'Pi{t) 






D,{t) 


to 




{T-(p} 


1 





t^ 




{T-ip} 





{r} 


to+T - (p 


F,n-l{T) 


{0} 


Fm{T - (fi) 


W} 


{to+T- <f) + 


F„_i(t) 





Fm{T - ip) 


{f} 


to+T 


f„(F„_i(T)+^) 





Fn-l{Fyn{T - if) 


{0} 


{t0+T) + 


F„(i?'„_i(r) + ^) 





Fr,-x{Fm{T - ^) + ^) 






Denote 4>i = Fn{Fm-i{r) + tp) and 02 = -F„_i(F^(r - (p) + ip). 

Proposition 2 shows that 0i = F„(F^_i(t) + </?)> F^+„_i(t) > F„_i(F^(t - (p) + (p) ^ 02- 
Then (^^(io + r + 1 - 0i) = 1, A(io + r + 1 - 0i) = 0, (/^^(io + r + 1 - ^i) = 02 + 1 - 0i and 
£)j(io + T + 1 — 0i) = 0. Thus, by Lemma 6, the pulse-coupled oscillator network can not achieve 
synchronization. □ 



5 Main theorem 

By using the lemmas above, we can obtain our main theorem. 

Theorem 1 Under the assumptions (Al) and (A2), from any initial phases {other than (pi{0) = 
■ ■ ■ — ip]s[{0)), the pulse-coupled oscillator network with delayed excitatory coupling can not 
achieve complete synchronization. 

Proof: If N — 2, then by Lemma 6 the theorem has been proved. Therefore, we only need to 
consider the case of A?" > 3. When A" > 3, for any initial phases (other than (pi{0) — ■ • • — (Pn{0)), 
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one of the following two cases must hold. 

Ceise 1: There exists a > such that the network is divided into two chques (A and B) 
at time ti. That is, the oscillators of clique A have been synchronized at ti and the oscillators 
of clique B have also been synchronized at ti, but the whole network has not been synchronized 
at ti. 

Because the initial phases do not satisfy (fii{0) — ■ ■ ■ — (fN{0), there must be a to with 
< to < ti such that the network is divided into two cliques (A and B) at time to, and the 
network is divided into three or more cliques at t' for all < t' < to- If to = 0, then the network 
is divided into two cliques at time 0. By Lemma 6, it can not achieve synchronization. If to > 0, 
similarly as in the proof of Lemma 5, there must be an oscillator i satisfying Di{to) = {0}. 
Then by Lemma 2(c), we have < (fi{to) < 1. Since the order of the oscillators is inessential 
to synchronization, we can re-arrange the order and assume that clique A includes oscillator 
1, . . . , m and clique B includes oscillator m + 1,. . . ,N, with 1 < m < A^. Without loss of 
generality, we assume that oscillator i is in clique A, i.e., 1 < i < m. By Lemma 3, we have 
< (pi{to) — ■ • ■ — iprn{to) < 1 and -Di(to ) — ■ ■ ■ — Dm{to) — 0. For clique B, one of the 
following two cases must hold. 

Subcase 1.1: There does not exist t" < to such that the oscillators of chque B have been 
synchronized at t". 

By Lemma 4, Di{to) = {0} shows -Dj(to) 7^ {0} for all m + 1 < j < A". Then by Lemma 5, 
we get 




(Pm{to) < 1 

• = DM = 0. 



So 




= CpNitt) = 



Similarly as in the proof of Lemma 6, the network can not achieve synchronization. 
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Subcase 1.2: There exists a t" < to such that the oscillators of clique B have been synchro- 
nized at t". 

By Lemma 4, we have Dm+i{to) = ■ ■ ■ — Disf{to)- There are two cases again. 

(i) Dm+l(to) ^ ■ ■ ■ ^ DN{to) ^ t 

Then we have 

Vi{to) = • • • = (Pm{to) < 1 
Di{t+) = ---^DUt+) = (l} 

V^m+l{to) = ■■■ = (pNito) < 1 

i^m+i(io) = --- = ^iv(io) = 
Denote (/?i(to) = • • • = </?m(^o) = 0i and <^m+i(^o) = • • • = </'iv(^o) = 02- H = ^2, then 

V'll^o) = ■ • • = iPmitQ) = 9?m+l(to) = • • • = 93 AT (to) = 01 < 1 

L'i(t+) = • • • - - Dm+i{4) = • • • = £'iv(4) = 0- 

By Lemma 3, the network has been synchronized at to, which contradicts that to < ti. Thus, 
01 7^ 02- Without loss of generality, we may assume 0i < 02- So, we have 

(pi{to + 1 - 02) = • • • = Vm{to + 1 - 02) = 01 + 1 - 02 < 1 

Di(to + 1 - 02) = • • • = D^(to + 1 - 02) = 

¥^^+1(^0 + 1 - 02) = • • • = (pN{to + 1 - 02) = 1 

L>^+i(to + 1 - 02) = • • • = D^ito + 1 - 02) = 
By Lemma 6, the network can not achieve synchronization. 

(ii) L'^+i(to) = • • • = D^ito) = {v} (A(to) = {0} shows < 77 < r). 
Then we have 

<Pl{to) = • • • = (Pm{to) < 1 

Di(t+) = ----Z^^(t+)=0 

Vm+l{to) = (^jv(to) < 1 

^m+l(to) = --- = ^iv(to) = M 

Denote (fii{to) — • • • — (fimito) — 0i and </?m+i(to) = • • • = (fiNito) — 02- It follows from 
Lemma 1 that 02 + 77 < 1. 
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li (f)i + Tj > 1, then 0i > 02- It follows that 

(pi{tQ + 1 - 0i) = • • • = (Pmito + 1 - 0l) = 1 

Di(to + 1 - 0i) = • • • = i?„(to + 1 - 0i) = 

Vm+l{to + 1 - 0l) = • • • = </?iv(^0 + 1 - 0l) = 02 + 1 - 01 < 1 

L'„+i(io + 1 - 01) = • • • = L>^(io + 1 - 0i) = {r? - (1 - 0i)} 
where < ?7 — (1 — 0i) < r. By Lemma 7, the network can not achieve synchronization. 
If 01 + 77 < 1, then 01 < 02- It follows that 

V^l{to + V) = ■■■ = ^m{to + V) = ^W-m(01 + ^7) 

(fim+lih + V) ^ ■ ■ ■ ^ ¥'N{to + V) ^ Fjv-m-l(02 + v) 
Dm+l{to + 7])^---^ DN{to +V) = {0}. 

Let X be a set of oscillators: oscillator k' E K ii and only if -Dfe'(io) — {0}- And let n be 
the number of the oscillators in K. It is obvious that clique B r\K — and chque A D K. 
Prom Proposition 1 and 0i,02 < 1, it follows that 0i = (Pi{to) — f~^{f{<^i{to)) + {n — l)e) 
and 02 = v5m+i(^o) = /"^(/(v'm+il^o )) +n€). Since A(to) = {0} and D„,+i{to) = {r]} (0 < 
T] < r), by Lemma 2(d), we have v9i(to ) > V^m+il^)- So, /(02) - /(0i) = [/(v'm+i(to )) + 
ne] — [/((^j(tQ )) + (n — l)e] < e. From the concavity assumption of function /, it follows that 
/(02 + v)- /(0i + ??)<£• Thus, Fiy_rn{(pi +'>])¥' -fiv-m-i(02 + v)- Then, similarly as in the 
proof of Lemma 6, the network can not achieve synchronization. 

Case 2: There does not exist h > such that the network is divided into two chques at time 

h. 

Assume that the pulse-coupled network can achieve synchronization. 

Then there must be a to > such that the network has been synchronized at to, but the 
network has not been synchronized at t for all < t < to- Similarly as in the proof of Lemma 5, 
there must be an oscillator i such that Di{to) = {0}. Because there does not exist ti > such 
that the network is divided into two cliques at time ti, there must be two other oscillators j and 
k such that 

oscillator i and j have not been synchronized at t for all t < to; 
oscillator j and k have not been synchronized at t for all t < to; 
oscillator k and i have not been synchronized at t for all t < to- 

21 



By Lemma 4, Di(to) = {0} shows that Dj(to) 7^ {0} and Dk(to) 7^ {0}. So, by Lemma 5, we 
have 

(Pj{to) = (pk{to) = 1 

Djito) = Dk{to) = 

However, ipi(to) = ^jito) = fkito) = 1 contradicts Di(to) = {0} by Lemma 2(c). 
So, in this case, the network can not achieve synchronization too. 

In summary, the theorem is proved. □ 

Theorem 1 shows that the presence of transmission delays can lead to desynchronization. 
However, one would say that the definition of complete synchronization in the theorem is awfully 
strong: it requires that all oscillators reach perfect synchronization in finite time; and say that 
usually synchronization only requires that the phase differences between oscillators converge 
to zero as time goes to infinity, i.e., asymptotical synchronization. In MiroUo and Strogatz's 
model this distinction is irrelevant, since if the phase differences are sufficiently small, perfect 
synchronization will be achieved after the next round of firings. But this is no longer the case 
when there exist transmission delays. Large numbers of numerical examples show that when 
there exist transmission delays, the phase differences can be small but there could still be big 
variations in the times to process the pulses in the pipeline. Simulation 1 in the following gives 
an example. Basing on the numerical result, we conjecture that asymptotical synchronization 
also cannot be achieved in pulse-coupled oscillator networks with delayed excitatory coupling. 

If our conjecture is true, then another problem arises: since neither complete synchronization 
nor asymptotical synchronization can be achieved, what does the model do? We numerically 
found that the model usually converges to two or more phased-locked clusters. It may be used 
to explain the phenomenon of clustering synchronization in some species of firefiies, which follow 
the phase advance model [2]. Simulation 2 in the following shows this. 

6 Numerical simulations 

In this section, we give two numerical simulations. 
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Consider the pulse-coupled oscillator network with delayed excitatory coupling: 

/(<^) = -/ exp{ln ^ + I (/ = 1.05) 
N = 100 
e = 0.001 

T = 0.1 

Obviously, it satisfies all the conditions of Theorem 1. In particular, we have /(2t) + Ne — 
0.5789 < 1. So, we conclude that the network can not achieve complete synchronization. 



Simulation 1 

In this simulation, we illustrate that even if all the initial phases are very close to each other, 
there could still be big variations in the times to process the pulses in the pipeline. An example 
with the initial phases chosen from a uniform distribution on (0, 0.01] is plotted in Fig.4. 
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Fig.4: Stroboscopic view on the phases ipi{tk) oi N = 100 oscillators plotted each time tk the first oscillator 
fires its fcth time, i.e., ipi{tk) = 1- The initial phases are chosen from a uniform distribution on (0, 0.01]. 
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Simulation 2 



In this simulation, wc investigate the clustering phenomena. The oscillators are initialized 
with a uniform random distribution on (0, 1]. Fig.5 shows that the network can converge to two 
or more phased-locked clusters. 
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Fig.5: Stroboscopic view on the phases (pi{tk) of iV = 100 oscillators plotted each time tk the first oscillator 
fires its fcth time, i.e., (pi{tk) = 1. The initial phases are chosen from a uniform distribution on (0, 1]. Eventually, 
the network is divided into two or more phased-locked clusters, for example, (a) two clusters (b) three clusters. 
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7 Conclusions 

In this paper, pulse-coupled oscillator networks with delayed excitatory coupling are studied. 
We propose an assumption (A2), which is reasonable to real biological systems, especially to 
fireflies. It is proved that under the assumptions (Al) and (A2), from any initial phases (other 
than 9?i(0) — ■ ■ ■ = ipn{0)), the network can not achieve complete synchronization. This result 
can explain why Photinus pyralis rarely synchronizes flashing, which is known as an example of 
pulsc-couplcd oscillator networks with delayed excitatory coupling. Furthermore, according to 
Simulation 1, we conjecture that asymptotical synchronization also cannot be achieved; and in 
Simulation 2, we exhibit a phenomenon of clustering synchronization. 
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